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ABSTRACT 


Quantitative  and  qualitative  predictions  of  galvanic  cor¬ 
rosion  rates  for  tube/tubesheet  galvanic  couples  are  presented 
and  discussed.  Quantitative  galvanic  corrosion  rate  predictions 
were  accomplished  using  a  long-term  electrochemical  potentiostatic 
polarization  data  and  finite  element  analysis.  The  120  day 
potentiostatic  polarization  data,  established  for  each  material 
in  natural  seawater,  were  utilized  to  incorporate  long  term 
polarization/time  effects.  Corrosion  rates  were  predicted  in 
terms  of  potential  and  galvanic  current  as  a  function  of  distance 
along  the  length  of  the  tube  from  the  tubesheet.  Qualitative 
galvanic  corrosion  predictions  were  made  by  examining  electro¬ 
chemical  polarization  data  and  geometry  in  the  context  of  the 
Wagner  polarization  parameter,  under  the  conditions  of  secondary 
and  tertiary  current  distribution. 

Seawater  exposures  of  each  tube/ tubesheet  arrrangement  are 
compared  and  contrasted  with  the  predicted  results.  These  ex¬ 
posures  used  segmented  tubes  connected  in  series  with  zero  re¬ 
sistance  ammeters.  Material  combinations  exposed  included 
a  Monel  alloy  400  tubesheet  coupled  to  a  90-10  copper-nickel 
tube,  a  90-10  copper-nickel  tube  coupled  to  a  zinc  anode  at 
the  position  of  the  tubesheet,  and  a  nickel-aluminum-bronze 
tubesheet  coupled  to  an  alloy  625  tube.  Analysis  was  also  con¬ 
ducted  on  a  nickel-aluminum-bronze  tubesheet  coupled  to  a  titanium 
50  tube,  and  comparisons  are  made  with  similar  data  reported  in 
the  literature. 

This  paper  is  directed  towards  demonstrating  the  need  for 
accurate  electrochemical  polarization  data  to  get  accurate  predic¬ 
tions.  The  feasibility  of  finite  element  analysis  for  accurate 
qualitative  as  well  as  quantitative  galvanic  corrosion  prediction 
for  the  tube/tubesheet  configuration  is  discussed. 
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INTRODUCTION 

Prediction  of  the  rate  of  galvanic  corrosion  by  the  use  of  galvanic  corrosion 
rate  tables  or  by  techniques  based  on  electrochemical  polarization  are  useful  but 


‘do  not  take  into  account  current  distribution  effects 


Thus,  these  methods  are  only 


IV 


Jn  the 


directly  applicable  when  there  simultaneously  exists^in  the  system  of  interest  1)  a 
1:1  anode  to  cathode  area  ratio,  2)  a  high  conductivity  electrolyte,  and  3)  a  sim¬ 
ple  geometry  with  exactly  known  electrochemically  active  surface  areas.  These  con¬ 
ditions  are  necessary  to  insure  that  the  current  is  uniformly  distributed  over  both 
the  anode  and  cathode.  For  most  actual  components  in  service  these  conditions  are 
not  satisfied,  and  simple  prediction  methods  such  as  those  above  cannot  be  applied. 


Ini  these  instances  the  distribution  of  galvanic  current  is  not  uniform. 


'A  common  geometric  configuration  that  exists  on  ships  where  dissimilar  metals 
are  used  and  where  galvanic  current  distribution  is  not  expected  to  be  uniform  is  the 
tube  and  tubesheet  configuration  in  a  condenser  or  heat  exchanger.  Higher  galvanic 
current  densities  are  expected  on  the  tube  surface  close  to  the  tubesheet,  with  the 
current  falling  to  insignificant  levels  at  some  distance  away  from  the  tubesheet. 

For  predicting  tube  corrosion  where  the  tube  is  more  active  than  the  tubesheet,  or 
for  predicting  tubesheet  corrosion  if  the  reverse  is  true,  the  distance  of  galvanic 
interaction  is  most  important  since  it  changes  the  effective  anode  to  cathode  ratio, 


which  has  a  significant  bearing  on  the  overall  corrosion  rate  of  the  anode. 
fty  r  o  t  V  6  f  rtktP  "Vd  .  {'St  ^ 


In  this  study  galvanic  corrosion  tendencies  of  tube/tubesheet  combinations  in 
seawater  were  analyzed  utilizing  long  term  potentiostatic  polarization  curves^ >2  as 
boundary  conditions  for  finite  elements  analysis  using  two  different  programs:  WECAN 
(Westinghouse  Electric  Computer  Analysis),  and  MARC.  Experimental  verification  was 
performed  using  segmented  tubes  electrically  coupled  to  tubesheets.  The  study  in¬ 
volved  several  condenser  tube/tubesheet  arrangements.  90-10  copper-nickel  tubing 
coupled  to  a  Monel  alloy  400  tubesheet,  alloy  625  tubing  coupled  to  a  nickel-aluminum- 
bronze  tubesheet,  and  a  90-10  copper-nickel  tube  coupled  to  anode  grade  zinc  at  the 
tubesheet  position  were  modelled  and  tested  under  quiescent  and  flowing  (2.4  m/s) 
conditions . 


In  addition,  predicted  galvanic  current  distributions  were  obtained  for  a 


titanium  (Ti-50)  tube  coupled  to  a  nickel-aluminum-bronze  tubesheet,  and  a  90-10 
copper-nickel  tube  coupled  to  a  simulated  Muntz  metal  tubesheet  represented  by  a 
fixed  potential  of  -250  mV  (versus  Ag/AgCl).  Seawater  exposures  for  comparison  in 
these  cases  were  obtained  from  available  literature.  All  bimetal  combinations  in¬ 
vestigated  are  summarized  in  Table  1. 


Table  1.  Tube/ tubesheet  finite  element  predictions. 


Tube 

Tubesheet 

Polarization 

Data 

Tube  Diameter 

Material 

Material 

Velocity  Time,  days 

cm 

Alloy  625 

Ni-Al-bronze 

2.4  m/s 

120 

2.66 

90-10  Cu-Ni 

Monel  alloy  400 

2.4  m/s 

120 

1.34 

Alloy  625 

Ni-Al-bronze 

2.4  ra/s 

1 

2.66 

90-10  Cu-Ni 

Monel  alloy  400 

2.4  m/s 

1 

1.34 

Alloy  625 

Ni-Al-bronze 

Quiescent 

120 

2.66 

90-10  Cu-Ni 

Monel  alloy  400 

Quiescent 

120 

1.34 

90-10  Cu-Ni 

Anode  grade 
zinc 

Quiescent 

30 

1.34 

Ti-50 

Ni-Al-bronze 

Quiescent 

1 

1.34 

Table  2.  Summary  of  tube/tubesheet  galvanic  corrosion  exposures 


Investigators 

Material  Combination 

Galvanic  Interaction 

Velocity 

Time 

Anode  Tubesheet 

Cathode  Tube 

Distance,  diameters 

m/  s 

Astley 

Ni-Al-bronze 

Titanium 

60 

1 

-600  mV  SCE 

70-30  Cu-Ni 

50 

quiescent 

1  hr 

18 

flowing 

2  hr 

35 

flowing 

13  d 

Gehring, 

-250  mV  SCE 

AL-6X 

120-480 

2.1 

Kuester , 

-250  mV  SCE 

AL-29-4 

120-480 

and  Maurer 

Muntz 

Ti-50A 

120-480 

Muntz 

90-10  Cu-Ni 

12-24  . 

Muntz 

Al-brass 

12-24 

Gehring 

Muntz 

Ti-50A 

120-240 

2.1 

Gehring 

Muntz 

A1-6X 

120-240 

2.1 

and  Kyle 

Al-bronze 

Ti-50A 

Gehring,  et  al. ,3~5  studied  a  "simulated"  Muntz  alloy  tubesheet  by  potentiostatic 
polarization  of  tubes  to  -250  mV  versus  SCE  at  the  tube  inlet.  Later,  actual  Muntz 
alloy  and  aluminum-bronze  alloy  tubesheets  were  coupled  to  more  noble  tube  materials 
in  flowing  seawater  at  7  ft/s  (2.1  m/s).  Noble  tube  materials  included  AL-6X  aus¬ 
tenitic  stainless  steel,  AL-29-4  ferritic  stainless  steel,  titanium  50  A,  90-10 
copper-nickel,  and  aluminum  brass.  For  both  simulated  and  actual  tube/tubesheet 
studies,  cathodic  polarization  of  tube  materials  along  the  length  of  the  tube  to  a 
distance  of  120  to  480  tube  diameters  was  observed  for  the  AL-6X,  AL-29-4,  and  Ti- 
50A  tubes.  For  couples  involving  both  90-10  copper-nickel  and  aluminum-brass 
coupled  to  a  muntz  metal  tubesheet,  cathodic  polarization  down  the  length  of  the 


tube  was  observed  for  a  distance  of  12  to  24  tube  diameters. 


D.J.  Astley  studied  a  titanium  tube  coupled  to  a  nickel-aluminum-bronze  plate 
in  seawater  flowing  at  1  m/s.^  The  length  of  titanium  tube  that  was  cathodically 
polarized  was  60  tube  diameters.  In  a  separate  study,  a  70-30  copper-nickel  tube  was 
potentiostatically  polarized  to  -600  mV  at  the  tube  inlet,  simulating  cathodic  pro¬ 
tection.  Under  quiescent  conditions  the  length  of  tube  receiving  cathodic  polariza¬ 
tion  after  one  hour  was  50  tube  diameters.  Under  flowing  conditions  after  two  hours, 
the  length  of  tube  receiving  cathodic  polarization  was  18  tube  diameters,  but  in¬ 
creased  to  35  diameters  after  13  days. 

In  summary,  the  galvanic  interaction  distance  from  all  studies  was  found  to  vary 
from  several  tube  diameters  to  several  hundred  depending  upon  the  materials'  polariza¬ 
tion  behavior,  the  geometry,  and  the  electrolyte. 

NUMBERICAL  ANALYSIS  FOR  GALVANIC  CORROSION 

Use  of  finite  element  and  boundary  element  techniques  has  shown  great  promise 
for  predicting  potential  and  current  distributions  for  galvanic  corrosion,  crevice 
corrosion,  and  cathodic  protection. 7-14  These  techniques  are  computationally  intense, 
and  require  powerful  computing  capabilities  for  all  but  the  most  simple  geometries. 
They  require  as  input  the  geometry  of  the  configuration  of  interest  and  the  con¬ 
ductivity  of  the  electrolyte  through  which  current  flows.  Boundary  conditions  on 
free  or  insulated  surfaces  are  that  current  cannot  flow  through  the  surface,  creating 
a  "reflection"  effect.  At  metal  surfaces,  the  polarization  characteristics  of  the 
metal  involved,  at  the  approprieate  flow  and  exposure  duration,  are  used  as  the 
boundary  conditions.  Output  from  these  analyses  are  potential  and  current  density  as 
a  function  of  position  along  the  structure. 

The  intent  of  this  study  is  to  demonstrate  that  tube/ tubesheet  galvanic  corrosion 
behavior  as  predicted  by  finite  element  analysis  depends  significantly  on  the  accuracy 
of  long  terra  polarization  data  and  open  circuit  potentials  used  for  each  material  in 
the  couple,  as  well  as  on  the  geometry  and  solution  conductivity. 


EXPERIMENTAL 


TUBE/ TUBESHEET  EXPOSURE  APPARATUS  AND  PROCEDURE 

Two  tube/tubesheet  combinations  were  tested,  both  at  quiescent  and  under  flowing 
(2.4  m/s)  conditions.  The  tubesheet  was  simulated  by  a  ring  of  material  mounted 
around,  but  insulated  from,  the  tube  end.  This  ring  was  coated  on  all  surfaces  ex¬ 
cept  that  surface  which  would  face  the  waterbox  in  a  real  condenser.  The  tube  was  cut 
into  segments,  each  insulated  from  the  others,  which  were  shortest  close  to  the  tube- 
sheet.  The  tube  segments  were  electrically  connected  to  each  other  and  to  the  tube- 
sheet  ring  with  zero-resistance  ammeters.  Holes  were  drilled  through  the  tube  wall 
periodically  to  facilitate  installation  of  micro-reference  electrodes  for  potential 
measurements.  Natural  seawater,  filtered  and  maintained  at  25°C,  was  placed  in 
contact  with  the  inside  of  the  tube  and  the  uncoated  face  of  the  tubesheet.  Cur¬ 
rent  and  potential  profiles  were  recorded  regularly  over  a  120  day  period.  At  the 
conclusion  of  the  quiescent  exposure  of  the  copper-nickel  tube,  this  tube  was  con¬ 
nected  to  a  sacrificial  zinc  anode  and  the  new  current  and  potential  profiles  mon¬ 
itored  over  an  additional  30  days  under  quiescent  conditions.  The  experimental  set¬ 
up  is  illustrated  in  Figs.  1  and  2.  The  two  material  combinations  used  were:  Alloy 
625  tube  with  Ni-Al -bronze  tubesheet,  and  90-10  Cu-Ni  tube  with  Monel  alloy  400 
tubesheet.  Dimensions  of  the  tubes  and  tubesheet  rings  are  listed  in  Table  3,  and 
the  dimensions  of  the  tube  segments  are  listed  in  Table  4. 

POLARIZATION  PROCEDURES  AND  REYNOLDS  NUMBER  COMPARISONS 

Data  from  120  day  polarization  experiments  on  condenser  materials  in  natural  sea¬ 
water,  as  described  elsewhere-*- >2  were  utilized  as  boundary  conditions  in  the  finite 
element  analysis  programs.  In  that  study,  individual  samples  were  held  at  a  constant 
potential  for  120  days.  Applied  currents  were  monitored  for  the  duration  of  the  ex¬ 
posure.  In  the  case  of  the  2.4  m/ s  flowing  seawater  polarization  data,  a  Reynolds 


number  of  30500  was  obtained,  indicative  of  turbulent  conditions.  For  the  1.34-cm- 
diameter-tube  configuration  described  above,  an  approximately  identical  Reynolds  num¬ 
ber  was  obtained.  For  the  case  of  the  2.66-cm-diameter  tube  described  above,  a  some¬ 
what  different  Reynolds  number  was  obtained,  however  turbulent  flow  was  still  main¬ 
tained.  In  addition  the  maximum  change  in  Nusselt  number  was  estimated  to  be  less 
than  1.75.  This  indicates  that  mass  transport  conditions  are  quite  similar. 

Table  3.  Dimensions  of  tube/tubesheet  condensers. 

In  625  Tube/Ni-Al-Bronze  Tubesheet  90  Cu-Ni  Tube/Monel  400  Tubesheet 
R  1.328  cm  (0.52  in.)  0.6693  cm  (0.2635  in.) 

Tl  0.0889  cm  (0.035  in.)  0.1245  cm  (0.049  in.) 

T2  0.5715  cm  (0.225  in.)*  0.2374  cm  (0.0937  in.)** 

Li  304.8  cm  (10  ft)  304.8  cm  (10  ft) 

L2  335.28  cm  (11  ft)  335.28  cm  (11  ft) 

*  Calculated  based  on  a  tube  center  to  tube  center  distance  of  3.9624  cm 
(1.56  in.)  on  a  triangular  pitch 

**  Calculated  based  on  a  tube  center  to  tube  center  distance  of  2.062  cm 
(0.812  in.)  on  a  triangular  pitch 

FINITE  ELEMENT  ANALYSIS 

Finite  element  analysis  was  conducted  by  John  Fu  at  Westinghouse  using  a  pro¬ 
gram  called  WECAN  on  both  material  combinations  tested,  but  for  quiescent  conditions 
only.  For  this  analysis  an  element  mesh  shown  in  Fig.  3  was  designed  using  axi- 
symraetrical  (cylindrical)  elements  to  reduce  the  three-dimensional  situation  to  a 
two-dimensional  analysis.  WECAN  models  polarization  by  a  series  of  straight  lines 
connecting  points  on  a  polarization  resistance  curve  derived  from  the  polarization 
curve  for  the  material  involved  at  the  appropriate  flow  rate.  The  polarization  re¬ 
sistance  is  the  overpotential  divided  by  the  current  density  at  that  overpotential, 
and  thus  looks  like  a  secant  line  on  the  polarization  curve  connecting  the  current 
at  the  freely-corroding  potential  to  the  current  at  the  potential  of  interest. 


Table  4.  Segmented  tube  sections  for  tube/tubesheet  galvanic 
corrosion  analysis. 

90-10  Copper  Nickel  Tube  Inconel  Alloy  625  Tube 

Coupled  to  Monel  Alloy  Coupled  to  Nickel-Aluminum 

400  Tubesheet  Bronze-Tubesheet 


Distance  from 
Tubesheet 
(cm) 

Tube  section 
Length 
(cm) 

Distance  from 
Tubesheet 
(cm) 

Tube  Section 
Length 
(cm) 

0 

1.0 

0 

3.0 

1.0 

2.0 

3.0 

3.0 

3.0 

3.0 

6.0 

3.0 

6.0 

5.0 

9.0 

3.0 

11.0 

5.0 

12.0 

4.0 

16.0 

25.0 

16.0 

4.8 

41.0 

50.0 

20.8 

5.0 

91.0 

63.8 

25.8 

5.0 

154.8 

75.0 

30.8 

5.0 

229.8 

75.0 

35.8 

8.0 

43.8 

8.0 

51.8 

8.0 

59.8 

9.0 

68.8 

10.0 

78.8 

10.0 

88.8 

10.0 

98.8 

10.0 

108.8 

26.0 

134.8 

30.0 

164.8 

30.0 

194.8 

40.0 

234.8 

70.0 

Finite  element  analysis  was  also  conducted  by  Raymond  Munn  at  the  Naval  Under¬ 
water  Systems  Center  using  a  program  called  MARC.  This  analysis  used  a  slightly 
different  axi -symmetric  mesh,  but  with  roughly  the  same  number  of  elements.  For  the 
MARC  analysis,  polarization  behavior  is  input  by  fitting  an  analytical  expression  to 
the  actual  polarization  data.  The  program  then  uses  the  potential  at  the  nearest 
integration  point  to  substitute  into  the  expression  to  calculate  the  current  density 
at  any  given  point  on  the  surface.  This  current  density  is  then  used  as  the  boundary 
condition  for  the  next  iteration  in  the  solution  of  the  problem. 


RESULTS 


ALLOY  625-Ni-Al-BRONZE  CONFIGURATION 

The  WECAN  analysis  for  the  120  day  quiescent  case  is  compared  to  120  day  ex¬ 
perimental  data  in  Fig.  4.  There  is  considerable  scatter  in  the  current  densities 
measured,  probably  due  to  crevice  corrosion  of  the  alloy  625  observed  at  the  con¬ 
clusion  of  the  test.  Considering  the  scatter,  agreement  is  good  between  predicted 
and  measured  current  densities.  Qualitative  agreement  of  potential  data  was  con¬ 
sidered  to  be  good  for  the  following  reason.  The  interaction  distance,  the  distance 
after  which  the  effect  of  the  tubesheet  on  potential  and  current  is  minimal,  is  sim¬ 
ilar  between  the  predicted  and  measured  data.  This  distance  is  about  150  cm,  which 
is  about  60-tube-diameters .  Therefore,  it  can  be  concluded  that  a  large  area  of 
cathode  is  present  in  this  case.  A  large  corrosion  rate  was  predicted  for  the 
Ni-Al-bronze  and  this  finding  was  actually  observed  once  the  cell  was  disassembled. 
There  is,  however,  some  difference  between  the  actual  potentials  observed  and  those 
predicted.  This  Is  because  the  polarization  data  utilized  in  the  analysis  does  not 
have  the  same  open  circuit  potential  for  the  alloy  625  as  the  tube  material,  thus 
shifting  the  whole  predicted  curve  below  the  measured  data.  Figure  5  illustrates 

the  potential  versus  time  behavior  for  alloy  625  in  quiescent  seawater.  Note  the 
fluctuations  possible  even  for  one  sample  in  quiescent  seawater.  These  fluctuations 
were  considered  to  be  an  indication  of  localized  corrosion  which  was  later  observed. 

The  MARC  analysis  for  the  120  day  quiescent  case  is  shown  in  Fig.  6.  Since  the 
same  polarization  data  was  used  in  this  analysis,  the  same  difficulty  in  the  predic¬ 
tion  of  potential  exists.  In  spite  of  this,  the  interaction  distance  is  about  60 
diameters,  in  agreement  with  the  WECAN  analysis  and  measured  potential  values  as 
well  as  the  measured  experimental  currents. 


Only  MARC  was  used  to  analyze  this  material  combination  under  flowing  (2.4  m/s) 
conditions.  Analyses  were  run  using  polarization  data  after  one  day  as  well  as  120 

days.  The  1  day  data  is  shown  in  Fig.  7  and  the  120  day  as  Fig.  8.  There  is  marginal  . 

qualitative  agreement  as  to  the  distances  of  galvanic  interaction  in  both  cases.  As 
with  the  quiescent  data,  the  freely-corroding  potential  of  alloy  625  used  in  the  analy¬ 
sis  appears  not  to  match  that  actually  observed.  Figure  8  shows  that  the  interaction 
distance  has  decreased  20  diameters  due  to  the  flow.  This  distance  is  similar  between 
the  computer  prediction  and  the  measured  potentials.  Although  the  1  day  prediction  in 
Fig.  7  (interaction  distances  of  15  diameters)  is  similar  to  the  120  day  prediction  in 
Fig.  8,  (interaction  distance  of  25  diameters)  the  measured  potentials  are  somewhat 
different.  In  Fig.  9  the  measured  tube/tubesheet  potential  profiles  can  be  seen  as 
they  develop  over  time.  The  galvanic  interaction  is  initially  low  (1  day)  when  the 
potential  difference  between  the  Ni-Al-bronze  and  the  alloy  625  is  small,  as  shown  in 
Fig.  10,  and  accordingly  the  galvanic  current  is  low.  When  the  open  circuit  potential 
of  the  alloy  625  increases,  the  interaction  distance  increases.  However,  the  galvanic  . 

corrosion  is  more  severe  at  70  days  in  this  case  than  with  the  quiescent  case  since 

the  cathodic  mass  transport  limiting  current  density  on  the  alloy  625  is  greater 
under  flowing  conditions  than  in  the  quiescent  case.  In  summary,  the  tube  potentials 
start  out  low  after  1  day,  then  increase  to  a  maximum  value  after  40  days.  Then  they 
decrease  until  70  days,  after  which  they  behave  erratically.  The  systematic  shifts 
in  tube/tubesheettref lect  the  time/potential  history  of  the  alloys  when 
freely  corroding,  as 

shown  in  Fig.  10.  After  1  day  the  freely  corroding  potentials  of  the  alloy  625  and 
the  bronze  are  similar  and  relatively  low,  but  the  alloy  625  potential  rises  after  a 
few  days,  creating  higher  galvanic  currents  and  higher  coupled  potentials.  After  day 
70  the  alloy  625  potential  may  have  behaved  erratically,  and  the  Ni-Al-bronze  has 


corroded  sufficiently  to  cause  the  remaining  copper  surface  to  exhibit  greater  pas¬ 
sivity,  as  observed. ^  This  explains  the  confusing  time  behavior  of  the  2.4  m/s  tube/ 
tubesheet  couple  reasonable  well  and  points  to  the  necessity  for  careful  consideration 
of  time  effects. 


90-10  Cu-Ni/MONEL  ALLOY  400  CONFIGURATION 

The  WECAN  analysis  of  this  material  combination  under  quiescent  conditions  at  120 
days  is  presented  in  Fig.  11.  The  analysis  of  this  material  combination  is  complicated 
by  the  similarity  in  freely-corroding  potentials  of  the  two  materials.  The  predicted 
interaction  distance  based  on  current  density  is  about  10  cm  or  8  diameters.  The 
measured  interaction  distance  was  closer  to  130  cm  or  100  diameters  after  120  days. 

The  explanation  for  this  behavior  seems  to  be  associated  with  the  gradual  passivation 
of  the  Cu-Ni  material  which  is  initially  anodically  polarized  near  the  tube  entrance. 
The  development  of  a  film  on  this  material  then  promotes  the  throwing  of  galvanic 
current  a  greater  distance  down  the  length  of  tube.  This  large  discrepancy  should  be 
considered  in  light  of  the  fact  that  all  of  the  current  is  at  least  a  factor  of  3  be¬ 
low  the  average  general  corrosion  current  for  freely-corroding  90-10  Cu-Ni,  and 
therefore  both  the  measured  and  predicted  galvanic  effects  are  insignificant.  The 
difference  between  predicted  and  observed  potentials  once  again  appears  large,  but 
all  potentials  lie  within  the  observed  range  of  freely-corroding  potentials  fot  the 
copper-nickel.  Thus,  the  differences  are  insignificant.  Since  the  increases  in  anode 
surface  area  minimizes  the  extent  of  galvanic  corrosion,  WECAN  was  successful  in  pre¬ 
dicting  an  insignificant  galvanic  corrosion  problem. 
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The  MARC  analysis  for  quiescent  120  day  conditions  in  Fig.  12  predicts  potentials 
similar  to  the  WECAN  analysis,  but  an  interaction  distance  more  in  agreement  with  that 
actually  measured.  The  reason  for  the  difference  between  the  two  analyses  is  not 
known,  but  MARC  was  successful  in  predicting  the  interaction  distance. 

Figures  13  and  14  present  the  MARC  analyses  and  measured  potentials  under  flowing 
conditions  for  1  day  and  120  days.  Potentials  were  different  between  predictions  and 
experimental  measurements  at  both  times.  In  addition,  predicted  glavanic  interaction 
distances  were  around  20  diameters  after  both  1  and  120  days,  whereas  the  observed 
distance  after  1  day  was  negligible  and  at  120  days  was  close  to  100  diameters.  As 
in  the  quiescent  case,  despite  these  differences,  an  insignificant  amount  of  galvanic 
corrosion  between  these  alloys  was  predicted  as  might  be  expected  from  the  small  open 
circuit  potential  differences  and  the  area  ratio  favorable  to  the  anode  material. 

Figure  15  illustrates  the  measured  tube  potential  profiles  as  they 
develop  over  time.  The  potentials  start  out  electronegative  at  1  day,  then  in¬ 
creases  so  that,  with  exception  of  the  60  day  data,  all  potentials  are  electropositive 
of  the  initial  values.  The  reason  for  this  can  be  seen  in  Fig.  10.  The  freely  cor¬ 
roding  potentials  for  both  alloys  increase  over  the  first  few  days  of  test.  The 
small  potential  difference  at  the  beginning  few  days  of  the  test  leads  to  an  insig¬ 
nificant  galvanic  couple  as  observed  and  predicted  for  the  1  day  case.  This  was 
reflected  in  the  measured  tube  potials  after  1  day,  where  no  interaction  distance 
could  be  observed.  MARC  was  successful  in  predicting  a  lack  of  significant  galvanic 
effect  under  these  circumstances.  Again,  exposure  time  effects  are  shown  to  be 
quite  substantial.  In  this  case  of  an  anodic  tube,  all  errors  were  on  the  conserva¬ 
tive  side,  indicated  by  a  large  anodecathode  area  ratio. 


90-10  Cu-Ni/ZINC  CONFIGURATION 


Figure  16  illustrates  the  results  of  the  MARC  analysis  for  90-10  Cu-Ni  tubing 
coupled  to  sacrificial  zinc  anode  material.  This  prediction  is  compared  to  the 
measured  values  from  this  study,  keeping  in  mind  that  the  measurements  were  on 
copper-nickel  which  had  already  been  exposed  for  120  days,  while  the  analysis  was 
based  on  polarization  data  from  only  a  30  day  exposure.  This  figure  also  contains 
measured  data  from  Astley^  on  a  1.1-cm  70-30  Cu-Ni  tube  cathodically  protected  at 
the  end  to  -600  mV.  Interaction  distances  were  50  diameters  fot  the  prediction,  50 
for  the  measurements  from  this  study,  and  40  for  Astley's  work.  This  agreement  is 
excellent,  considering  the  differences  in  the  various  studies. 

OTHER  MATERIAL  COMBINATIONS 

Figure  17  illustrates  a  potential  profile  predicted  for  a  1 .34-cm-diameter 
titanium  tube  coupled  to  a  Ni-Al-bronze  tubesheet.  The  MARC  prediction  is  based  on 
quiescent  polarization  behavior  after  30  days,  and  the  experimental  measurements  are 
for  D.J.  Astley's  1 .46-cm-diameter-tube  with  1  m/ s  flow. 6  Predicted  potentials  are 
quite  different  from  those  measured.  This  discrepancy  Is  probably  due  to  a  large 
change  in  the  polarization  characteristics  of  titanium  when  going  from  quiescent 
conditions  to  a  flow  of  1  m/s.  The  prediction  is  for  an  interaction  distance  of  over 
100  diameters  whereas  the  measured  distance  was  60  diameters,  which  is  reasonable 
agreement.  In  both  cases  a  large  cathode  to  anode  area  ratio  is  predicted  and 
galvanic  corrosion  of  the  bronze  is  expected  to  be  severe. 

Figure  18  shows  a  predicted  potential  distribution  on  a  134-cm  diameter  90-10 
Cu-Ni  tube  at  2.1  m/ s  coupled  to  a  constant  potential  source  of  -250  mV,  simulating 
a  Muntz  metal  tubesheet.  The  polarization  data  used  for  the  potential  predictions 
was  obtained  at  2.4  ra/ s  flow  at  1  day.  The  measured  current  distribution  data  is 


from  work  by  Gehring,3~5  who  did  not  report  potential  distribution.  The  predicted 
interaction  distance  and  the  measured  distance  are  essentially  identical  at  about 
15  diameters.  In  this  case  galvanic  corrosion  is  limited  to  the  tube  entrance 
despite  the  presence  of  the  large  cathodic  tube  surface. 

DISCUSSION 

QUANTITATIVE  PREDICTION  BY  THE  USE  OF  FINITE  ELEMENT  ANALYSIS 

As  described  in  the  preceeding  section,  agreement  between  predicted  and  measured 
potential  profiles  is  generally  not  very  good.  This  may  be  due  to  sensitivity  in  the 
computer  analysis  to  small  changes  in  the  shape  of  the  polarization  curves  used  for 
boundary  conditions.  Also,  the  discrepancy  may  be  due  to  the  variable  nature  of  the 
freely-corroding  potential  of  a  material  in  seawater,  which  can  lead  to  a  displace¬ 
ment  of  the  potential  profiles  but  seems  to  affect  the  shape  of  these  profiles  only 
slightly. 

Table  5  presents  a  summary  of  predicted  versus  measured  interaction  distances 
for  all  tube/tubesheet  combinations  discussed  above.  Interaction  distances  were 
measured  to  range  from  15  to  100  diameters,  and  predicted  to  range  from  8  to  100 
diameters.  With  the  exception  of  one  WECAN  analysis,  measured  and  predicted  inter 
action  distances  were  very  close  under  quiescent  conditions.  With  the  exception  of 
one  MARC  analysis,  measured  and  predicted  interaction  distances  under  flowing  con¬ 
ditions  were  also  quite  similar.  Both  exceptions  occurred  in  the  tube/tubesheet 
combination  of  90-10  Cu-Ni/  Monel  alloy  400,  where  the  total  amount  of  galvanic 
interaction  is  extremely  small  due  to  similarity  in  the  materials'  freely-corroding 
potentials  and  the  passive  behavior  of  these  alloys  in  seawater.  Thus,  the  finite 
element  technique  is  a  good  way  to  predict  the  interaction  distance  down  tubes,  and 


thereby  the  effective  area  ratio  of  a  tube/ tubesheet  galvanic  couple.  This  was  the 
case  over  a  range  of  material  combinations,  exposure  periods,  flow  velocities,  com¬ 
puter  programs,  and  data  sources. 

There  is  insufficient  information  in  this  study  to  make  any  statement  about 
prediction  accuracy  for  quantitative  determination  of  glavanic  corrosion  rates. 

Table  5.  Distances  of  maximum  galvanic  interaction. 


Tube 

Tubesheet 

Technique 

Interaction  Distance 

(Diameters ) 
Quiescent  2.4  m/ s 

Time 

(Days) 

Alloy  625 

Ni-Al-bronze 

Measured 

60 

15 

120 

25 

1 

WE  CAN 

60 

120 

MARC 

60 

15 

120 

15 

1 

90-10  Cu-Ni 

Monel 

Measured 

100 

100 

120 

WE  CAN 

8 

120 

MARC 

100 

20 

120 

20 

1 

90-10  Cu-Ni 

Zinc 

Meas  (Scully) 

50 

30 

Meas  (Astley) 

40 

1  h 

WE  CAN 

50 

30 

Ti-50 

Ni-Al-bronze 

Meas  (Astley) 

60 

WECAN 

100 

30 

90-10  Cu-Ni 

-250  mV 

Meas  (Gehring) 

15  (2.1  m/s) 

WECAN 

15 

1 

QUALITATIVE  PREDICTION  BY  THE  USE  OF  WAGNER  NUMBER  ANALYSIS 

The  Wagner  number^-5,16  can  be  used  to  evaluate  the  relative  degree  of  uniformity 
of  current  distribution  between  anodes  and  cathodes  such  as  galvanic  current.  This 
parameter  describes  the  ratio  of  the  kinetic  resistance  (polarization)  to  the  ohmic 
resistance  (ionic  conduction)  throughout  the  electrolyte  separating  the  anode  and  the 


Kinetic  Resistance 
Ohmic  Resistance 


cathode  as  follows: 


When  Che  Wagner  number  gets  very  large  a  uniform  current  distribution  over  the 
anode  and  cathode  is  promoted.  When  the  Wagner  number  gets  very  small  a  nonuniform 
current  distribution  is  promoted.  The  Wagner  number  can  be  changed  by  modifying 
either  the  polarization  characteristics  of  the  anode  and  cathode  materials  involved, 
or  the  solution  conductivity,  or  both,  thus  changing  the  uniformity  of  current  dis¬ 
tribution.  The  Wagner  number  analysis  is  applicable  to  primary,  secondary,  or 
tertiary  current  distribution.  For  the  primary  case,  only  the  solution  conductivity 
and  a  characteristic  length  between  anode  and  cathode  are  considered: 

W  »  K/L 

where  K  is  the  solution  conductivity,  and  L  is  a  characteristic  length.  In  this 
study  the  Wagner  number  was  utilized  under  conditions  of  secondary  and  tertiary  cur¬ 
rent  distribution,  meaning  that  linear  Tafel,  and  concentration  type  polarization 
conditions  were  all  considered  as  below: 

W  =  K  (5  rj/ 5i) lavg/L 

Table  6  illustrates  these  various  conditions.  Note  that  in  the  secondary  cur¬ 
rent  distribution  case  a  large  Tafel  coefficient  B  or  B'  promotes  a  large  Wagner  num¬ 
ber,  and  uniformity  of  current  distribution.  In  the  case  of  concentration  polariza¬ 
tion,  currents  near  the  limiting  current  density  promote  large  Wagner  numbers. 

Figures  19  and  20  summarize  the  various  cases  of  large  and  small  Wagner  numbers  on 
the  galvanic  current  distribution  for  both  the  case  of  an  anode  and  a  cathode  as  the 
tube  material,  respectively.  The  exact  magnitude  of  the  galvanic  current  is  not 
determined  from  the  type  of  treatment.  The  quantitative  solution  to  these  types  of 
problems  requires  the  use  of  numerical  methods  as  discussed  above. 
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Table  6.  Wagner  polarization  parameter  determination. 
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Relationship 
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Primary 

— 

— 

K/L 

Secondary 

n  =Bi 

( linear) 

B 

KB/L 

Secondary 

r,  =B’ln(i/io) 

(Tafel) 

B’  /i 

KB’/iL 

Tertiary 

rj  =  -RT/nF  (ln(l-i/iL) 

(diff .) 

RT/nF  l/(iL-i) 

KRT/nF(iL-i)L 

Figure  21  is  the  cathode  polarization  curve  for  alloy  625  in  quiescent  seawater. 
The  anodic  curve  for  nickel-alurainum-bronze  intersects  this  curve  in  the  diffusion 
limited  region,  which  would  indicate  that  the  Wagner  number  should  be  large  for  the 
tube  surface,  approaching  infinity.  A  large  Wagner  number  should  indicate  current 
uniformity,  which  should  result  in  a  large  interaction  distance  on  the  alloy  625. 
Observed  interaction  distances  are  in  fact  quite  large,  60  diameters.  Since  titanium 
has  similar  polarization  characteristics  to  alloy  625,  similar  behavior  should  be 
observed  with  couples  using  titanium  tubes.  As  Table  5  illustrates,  Astley  measured 
an  equivalent  60diaraeter  interaction  distance.  Under  flowing  conditions,  the  dif¬ 
fusion  limited  portion  of  the  cathodic  curve  for  alloy  625  should  be  limited  to  more 
more  electronegative  potentials  causing  the  Wagner  number,  and  the  resultant  inter¬ 
action  distance,  to  shrink.  This  was  actually  observed,  as  the  interaction  distance 
was  only  15-25  diameters  at  2.4  m/s  flow.  Thus,  the  Wagner  number  gives  a  good  qual¬ 
itative  indication  in  this  case  of  the  amount  of  tube  surface  that  has  significant 
effect  in  the  galvanic  couple. 
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Figure  22  is  the  cathodic  polarization  and  Wagner  analysis  of  copper-nickel  in 
quiescent  natural  seawater.  A  galvanic  couple  between  a  zinc  anode  at  the  position 
of  the  tubesheet  and  a  Cu-Ni  tube  would  behave  similarly  to  the  above  example  in 
quiescent  water  and  produce  a  large  Wagner  number,  since  cathodic  currents  on  the 
tube  are  near  the  limiting  current  density  for  oxygen  reduction.  Large  interaction 
distances  are  predicted  by  large  Wagner  numbers,  and  were  observed  at  40-50  diameters 
both  in  this  study  and  by  Astley,  as  seen  in  Table  5. 

A  Cu-Ni  tube  coupled  to  a  Muntz  metal  or  similar  brass  tubesheet  would  result 
in  a  low  Wagner  number  since  the  cathodic  current  on  the  Cu-Ni  is  far  below  the 
limiting  current  density.  This  is  seen  in  Fig.  22,  where  the  Wagner  number  is  rela¬ 
tively  low  at  potentials  of  -250  mV  or  above.  Table  5  shows  that  Gehring  reported 
only  a  15  diameter  interaction  distance  under  these  circumstances.  Again,  the  Wagner 
number  analysis  is  a  good  one  for  this  situation. 

The  couple  between  a  Cu-Ni  tube  and  a  Monel  alloy  400  tubesheet  will  result  in  a 
low  Wagner  number,  as  illustrated  schematically  in  ligs.  22  through  24.  This  should 
result  in  galvanic  current  being  limited  to  the  entrance  of  the  anodic  tube.  Low 
interaction  distances  are  predicted  in  this  case,  but  in  this  study  the  observed  in¬ 
teraction  distances  are  very  large,  100  diameters.  This  can  only  be  explained  if  the 
tube  material  becomes  passivated,  in  which  case  the  Wagner  parameter  becomes  large 
for  the  anode  tube  material  and  the  current  becomes  more  uniform. 

The  experimental  results  not  modelled  here,  but  reporte '  in  Table  2  can  also  be 
analyzed  using  the  Wagner  number.  Following  the  same  logic  used  above,  a  large  in¬ 
teraction  distance  is  predicted  for  all  studies  in  this  table  except  those  involving 
copper-nickel  or  brass.  With  only  2  exceptions  from  Astley's  work,  the  interaction 
distances  were  all  observed  to  be  from  50-480  diameters,  in  agreement  with  the  pre¬ 
diction.  The  copper-nickel  and  brass  couples  are  all  predicted  to  have  short 

interaction  distances,  and  all  were  reported  in  the  range  of  12-24  diameters. 
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The  above  analyses  indicate  how  good  a  predictor  the  Wagner  number  is  for  esti¬ 
mating  in  teraction  distances  in  tube/tubesheet  galvanic  couples.  Although  quan¬ 
titative  prediction  is  not  possible  using  this  analysis,  qualitative  information  can 
be  gained,  which  is  useful  for  determining  the  effects  of  exposure  variables  on  in¬ 
teraction  distance.  This  can  be  accomplished  by  examining  the  effects  of  these 
variables  on  polarization  behavior  and  determining  Wagner  numbers  as  a  function  of 
potential . 


CONCLUSIONS 

Quantitative  predictions  of  galvanic  corrosion  by  finite  element  methods  can  be 
accomplished  by  predicting  distances  of  galvanic  interaction  provided  that  long  term 
electrochemical  polarization  data  and  open  circuit  potential'  behavior  are  carefully 
considered  and  incorporated  into  such  analysis.  Without  consideration  of  potential- 
current-time  effects  resulting  predictions  may  be  misleading  or  in  error.  Qualitative 
indications  of  galvanic  corrosion  interaction  distances  may  be  obtained  through  utili¬ 
zation  of  the  Wagner  number  analyzed  in  the  case  of  the  primary  or  secondary  current 
distribution.  Consideration  of  potential-current-time  effects  are  still  required  and 
no  quantitative  data  may  be  yielded  through  this  method. 
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Fig.  6.  Galvanic  corrosion  behavior  for  an  alloy  625  tube 
coupled  to  Ni-Al -bronze :  quiescent  120  days. 


Fig.  7.  Galvanic  corrosion  behavior  for  alloy  625  tube  coupled 
to  Ni-Al-bronze  tubesheet:  flowing  1  day. 
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Galvanic  corrosion  behavior  for  alloy  625  tube  coupled 
to  Ni-Al-bronze  tubesheet:  flowing  120  days. 
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Galvanic  corrosion  behavior  time  effects  for  alloy  625  tube 
coupled  to  Ni-Al-bronze  tubesheet:  flowing  1  to  120  days. 
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Fig.  10.  Potential  versus  time  behavior  for  tube  and  tubesheet 
materials  in  flowing  seawater. 
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Fig.  11.  Current  and  potential  distribution  on  90-10  Cu-Ni  tube 
after  120  days  in  quiescent  seawater. 
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Fig.  12.  Galvanic  corrosion  behavior  for  a  90-10  Cu-Ni  tube  coupled 
to  Monel  alloy  400:  quiescent  120  days. 
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Galvanic  corrosion  behavior  for  a  90-10  Cu-Ni  tube  coupled 
to  Monel  alloy  400:  flowing  1  day. 
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Fig.  14.  Galvanic  corrosion  behavior  for  a  90-10  Cu-Ni  tube  coupled 
to  Monel  alloy  400:  flowing  120  days. 
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Fig.  15.  Galvanic  corrosion  behavior  time  effects  for  90-10  Cu-Ni  tube 
coupled  to  Monel  alloy  400  tubesheet:  flowing  1  to  120  days. 
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Fig.  18.  Galvanic  corrosion  behavior  for  a  90-10  Cu-Ni  tube  coupled  to 
Muntz  metal  alloy  tubesheet. 
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Fig.  20.  Schematic  galvanic  current  versus  distance  plot  for  a  cathode  tube 
coupled  to  an  anode  tubesheet  for  various  magnitudes  of  the  Wagner 
polarization  parameter. 


Fig.  21.  Schematic  E  ve  sus  log  i  plot  describing  galvanic  corrosion 
behavior  of  alloy  625  coupled  to  Ni-Al-bronze, 
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